Superconducting Gap Anisotropy and Quasiparticle Interactions: a Doping 

Dependent ARPES Study 
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Comparing ARPES measurements on Bi2212 with penetration depth data, we show that a description 
of the nodal excitations of the d-wave superconducting state in terms of non-interacting quasiparticles is 
inadequate, and we estimate the magnitude and doping dependence of the Landau interaction parameter 
which renormalizes the linear T contribution to the superfluid density. Furthermore, although consistent 
with d-wave symmetry, the gap with underdoping cannot be fit by the simple cosfcz-cosfcy form, which 
suggests an increasing importance of long range interactions as the insulator is approached. 
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There is little doubt about the fundamental impor- 
tance of many-body interactions in high temperature 
cuprate superconductors Q. Quantifying these interac- 
tions is difficult in the normal state of these materials, 
given the lack of well-defined single-particle excitations 
as revealed by various experiments. On the other hand, 
well-defined quasiparticle excitations do exist in the su- 
perconducting state, and it is believed that a description 
of the low temperature state in terms of superfluid Fermi 
liquid theory is appropriate. In Fermi liquid theory, the 
quasiparticles are characterized by a renormalized Fermi 
velocity Vp, and their residual interactions described by 
Landau parameters, which manifest themselves through 
a renormalization of various response functions relative 
to that given by a non-interacting theory. For example, 
in the cuprates, the Fermi velocity vp has been deter- 
mined by angle resolved photoemission (ARPES) stud- 
ies in Bi 2 Sr 2 CaCu20 8+5 (Bi2212) @ to be renormalized 
by a factor of two to three over that given by band the- 
ory. The strong renormalization of the superfluid density 
/? s (0) has also been known for some time, where one sees 
a scaling with the number of doped holes: the Uemura 
relation ||. 

In this paper we examine an issue which is at the heart 
of the nature of quasiparticles in the superconducting 
state of the cuprates, that is, whether the slope of the su- 
perfluid density at low temperatures, dp s /dT, is affected 
by interactions or not, and what the relation of its renor- 
malization is to that of /9 S (0), questions of considerable 
debate in the recent literature The importance of 

p s (T) to an understanding of cuprate superconductivity 
derives from the early observation of a linear T suppre- 



sion of p s (T) H , since this is explained most naturally by 
the thermal excitations of quasiparticles near the nodes 
of a d-wave superconducting gap. Related to this is the 
interesting question of whether the gap around the node 
scales with T c , as has been suggested from a recent anal- 
ysis of magnetic penetration depth data 0. 

To address these issues we use the unique capability of 
ARPES to directly measure the Fermi wavevector kp , ve- 
locity Vp, and the superconducting gap anisotropy near 
the node, from which we can estimate the slope of p s (T) 
assuming non- interacting quasiparticles. Comparing this 
with the actual value obtained by penetration depth ex- 
periments leads to a direct estimate of the renormaliza- 
tion due to quasiparticle interactions. This is done by 
exploiting the relation 
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where A is the penetration depth, and A is a doping- 
independent constant: A = Ahi2akBn/cd with a the 
fine structure constant, ks the Boltzmann constant, c 
the speed of light, and n the number of CuOi layers 
(4 for Bi2212) per c-axis lattice constant d (30.9 A for 
Bi2212). ARPES is used to determine the three pa- 
rameters at the node: the Fermi velocity Vp, the Fermi 
wavevector kp, and the slope of the superconducting gap 
v a = l/2\dA/d<f>\(4> — 7r/4), where 4> is the Fermi surface 
angle. The latter is normalized such that va — A max for 
the simple d-wave gap A(0) = A max cos(2</>). 

The only unknown in Eq. 1 is the renormalization fac- 
tor (3 due to quasiparticle interactions; in the isotropic 
Fermi liquid theory (3=1 + F\ s /2, where Fu is the 



1 



/ = 1 spin symmetric Landau parameter, and quantifies 
the backflow of the medium around the quasiparticles || . 
By comparing ARPES and penetration depth data, we 
estimate f3 and its doping dependence. In particular, dif- 
ferent assumptions in the recent literature about 
the doping dependence of va has led to different conclu- 
sions regarding the value and doping dependence of (3 in 
Eq. 1. 

Our main results are as follows. (1) We determine the 
doping dependence of the gap anisotropy from ARPES. 
Although consistent with a node on the Fermi surface 
along the zone diagonal (0 = 7r/4) for all doping levels, 
the shape of the gap changes with underdoping: while its 
maximum value increases we find the new result 

that the gap becomes flatter near the nodes, i.e. v\ de- 
creases. (2) Using our data on the doping dependence of 
i>a , we exploit Eq. 1 and use available values of the pene- 
tration depth A(T) |p|,^2|-^5| to estimate the renormaliza- 
tion factor f3. We find that (3 is considerably smaller than 
unity and decreases with underdoping, in contrast to pre- 
vious suggestions in the literature (3) Our results 
on the doping dependence of the gap anisotropy and its 
relation to penetration depth data provide important ev- 
idence that the strength of both the pairing interaction 
and the quasiparticle interactions increase with reduced 
doping. 

The ARPES experiments were performed at the Syn- 
chrotron Radiation Center, Wisconsin, using both a high- 
resolution 4-meter normal incidence and plane grating 
monochromators, with a resolving power of 10 4 at 10 11 
photons/sec. We used 22 eV photons, with a 17 meV 
(FWHM) energy resolution, and a momentum window 
of radius 0. 045-77 (in units of 1/a where a is the Cu-Cu 
separation). The high quality single crystal samples were 
float-zone grown, with doping changed by varying the 
oxygen partial pressure during annealing. All samples 
show sharp x-ray diffraction rocking curves and flat sur- 
faces after cleaving as determined from specular laser re- 
flection. We label the samples by their doping (UD for 
underdoped, OD for overdoped) and onset T c . 

Fig. 1 shows ARPES data at T=15 K for an UD75K 
sample at different k-points along the Fermi surface, kp 
was carefully chosen using the criterion that the leading 
edge of the spectrum has minimum binding energy with 
the steepest slope, when compared with other spectra 
along a cut perpendicular to the Fermi surface, as dis- 
cussed earlier fl6f| . The zero of binding energy (Ep) was 
determined from the spectra (not shown) of a polycrys- 
tallinc Pt reference in electrical contact with the Bi22I2, 
recorded at regular intervals to ensure accurate determi- 
nation of the Fermi energy, Ep. From the shift of spectral 
weight away from Ep, one clearly sees an anisotropic gap, 
which is maximal near the (n, 0) point (cf> — 0) and zero 
near the (tt, tt) direction (</> = 45°) 
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FIG. 1. Spectra of an UD75K Bi2212 sample (solid line) 
in the vicinity of Ef taken at T=15K, each labeled by the 
Fermi surface angle <f). For two angles we also plot spectra 
from an OD87K sample (dotted line). 
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FIG. 2. Values of the superconducting gap as a function of 
the Fermi surface angle <f) obtained for a series of Bi2212 sam- 
ples with varying doping. Note two different UD75K samples 
were measured, and the UD83K sample has a larger doping 
due to aging ]lq |. The solid lines represent the best fit using 
the gap function: = A max [B cos(20) + (1 — B) cos(60)] 
as explained in the text. The dashed line in the panel of an 
UD75K sample represents the gap function with B=l. 
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For comparison we also plot (dashed line) in Fig. 1 
ARPES spectra from an OD87K sample at two points 
on the Fermi surface. (For more OD data see Ref. Jl8[.) 
We immediately see that the UD sample has a larger 
maximum gap (<j) — 0) than the OD one, but it has a 
smaller gap at the corresponding point ((f) — 38 degrees) 
near the node. Thus the raw data directly give evidence 
for an interesting change in gap anisotropy with doping. 

To quantitatively estimate the gap, we have modeled 
the low temperature data by a simple BCS spectral func- 
tion, taking into account the measured dispersion and the 
known energy and momentum resolutions. Details of this 
analysis, and error estimates, have been described earlier 
in the context of OD samples jl9|,[l8| . The resulting an- 
gular dependence of the gap is plotted in Fig. 2 for six 
samples. 
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FIG. 3. (a) Normalized slope of the gap at the node 
(t>A/A max ) vs gap maximum A max . Note the clear drop from 
unity as one enters the underdoped regime, (b) Slope of the 
superfluid density (in units of 10" 9 A" 2 A" _1 ) vs T™ ax -T c es- 
timated from ARPES measurements based on non-interacting 
quasiparticles in Bi2212 (filled circles) compared with direct 
penetration depth measurements in YBCO (open- J13] and 
filled- @ triangles, T™ ax =92K) and Bi2212 (open- |jj and 
filled- [14] squares, T™ ax =95K). The error bars for the latter 
values were ±15% based on ±5% error bars for A(0) J^j. The 
lines are guides to the eye. 



To further quantify this change in anisotropy, we have 
used the following expression to fit the gap: Ak = 
A max [ J Bcos(20) + (l-S)cos(60)] with < B < 1, where 
B is determined for each data set. Note that cos(60) 
is the next harmonic consistent with ci-wave symmetry. 
We find that while the overdoped data sets are consistent 
with B « 1, the parameter B decreases significantly in 
the underdoped regime. To emphasize the significance of 
B < 1, we plot in the panel of an UD75K sample of Fig. 2 
a dashed curve with B = 1 along with the best fit curve 
for that sample. From these fits, one easily determines 



the value «a discussed earlier in the context of Eq. 1. 
In Fig. 3a, we plot WA/A max for seven samples (the six 
analysed above plus an UD85K sample from Ref. po|). 
One can clearly see from this figure the trend that under- 
doping leads to an increase in the maximum gap together 
with a decrease in the gap slope at the node. 

Several questions need to be addressed before proceed- 
ing further. First, could the flattening at the node be, 
in fact, evidence for a "Fermi arc" (a line of gapless ex- 
citations), especially since such arcs are seen above T c in 
the underdoped materials [pp|? Given the error bars on 
gap estimates in Fig. 2, it is impossible to rule out arcs 
in all the samples. Nevertheless, it is clear that there are 
samples (especially OD87K, UD80K and UD75K) where 
there is clear evidence in favor of a point node rather 
than an arc at low temperatures. Furthermore, it is very 
important to note that a linear T dependence of p s (T) at 
low temperature, for all doping levels, in clean samples 
gives independent evidence for point nodes |]6|, |l2| , ^5| . 

Second, is the change in gap anisotropy intrinsic, or 
related to impurity scattering? We can eliminate the lat- 
ter explanation on two grounds. The maximum gap in- 
creases as the doping is reduced, opposite to what would 
be expected from pair breaking due to impurities. Also, 
impurity scattering is expected to lead to a characteris- 
tic "tail" to the leading edge pl| , for which there is no 
evidence in the observed spectra (see Fig. 1). 

We suggest that the change in the gap function with 
underdoping is related to an increase in the range of the 
pairing interaction: the cos(60) term in the Fermi sur- 
face harmonics can be shown to be closely related to 
the tight binding function cos(2/c2;) — cos(2k y ), which 
represents next nearest neighbours interaction, just as 
cos(20) is closely related to the near neighbor interac- 
tion cos(fc K ) — cos(ky). On very general grounds, the 
increasing importance of the cos(6</)) term with under- 
doping could arise from a decrease in screening as one 
approaches the insulator. Similar effects also arise in 
specific models. In models of spin-fluctuation mediated 
d-wave pairing, an increase in the antiferromagnetic cor- 
relation length with underdoping leads to a more sharply 
peaked pairing interaction in k-space, causing a flatten- 
ing of the gap around the node as we find here. In inter- 
layer tunneling models, one also expects changes in the 
shape of the gap which might be correlated with doping 

We note that the ratio of the dispersion normal to the 
Fermi surface (vf) to that along the Fermi surface (v a) 
is quite large, 20 in the overdoped case, and becomes 
even larger as the doping decreases, in contrast to the 
undoped insulator which exhibits an isotropic dispersion 
about the (7r/2,7r/2) points p3| . This implies that the 
electronic dispersion in the superconductor in this region 
of the zone may not be as closely related to the insulator 
as has been recently suggested p4j . 

We now return to Eq. 1. It is known from previous 
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ARPES measurements that the band dispersion along 
(0,0) — (tt, 7t) is rather strong and doping independent 
p5| with an estimated vf — 2.5 x 10 7 cm/sec [jl9|. It 
is also known that kp along this direction is 0.737 AT 1 
and relatively doping independent |l6{ ]. Using these in- 
puts, together with the strongly doping dependent va, 
we can estimate the slope |e?A~ 2 /dT| in the case of non- 
interacting quasiparticles (/3 = 1). As shown in Fig. 3b 
(filled circles) we find the resulting slope is of order 
6 x 10~ A~ 2 K _1 and is reduced by approximately 30% 
in going from UD75K to OD87K. 

Fig. 3b also shows the values of |<iA~ 2 /e?T| obtained 
from London penetration depth measurements |^|,|l2] |l5| . 
Although, there is considerable variation in the mea- 
sured values of A(0) and low temperature dX/dT from 
one group to another, probably due to the use of different 
techniques, we find evidence for the following trend: the 
slope dp s /dT decreases with underdoping. For YBCO 
this effect is weak in the UBC data fl^| , but much 
stronger in the Cambridge data [fl5|| . The limited data 
available for Bi2212 are consistent with this trend |fl4) . 
The striking feature is that, in all cases, this trend in 
dp s /dT is exactly the opposite of that deduced from a 
theory with non-interacting quasiparticles {(3 = 1) using 
ARPES input. That is, from Fig. 3b, it is clear that 
the renormalization factor (3 is considerably smaller than 
unity and doping dependent, a conclusion different from 
that inferred earlier To get an estimate of the 

doping dependence of (3, we use the Bi2212 values of 
Ref. [jl4| for OD85K and UD80K samples in comparison 
to our own values on OD87K and UD80K, obtaining a (3 2 
of 0.32 and 0.17, respectively. This is roughly consistent 
with a [3 which varies as x, the number of doped holes, 
which would be the expected result from the x scaling 
of p s (0) |§JJ§. On the other hand, as noted earlier |^|, a 
weaker doping dependence of (3 seems to be implied by 
the UBC data. Given the difficulties of measuring the 
superconducting gap in YBCO by ARPES, this points 
to the need for further penetration depth experiments on 
Bi2212 samples so that a more detailed comparison to 
ARPES data can be made. 

In conclusion, we find that the gap anisotropy of 
Bi2212 changes strongly as a function of doping, implying 
an increase in the range of the pairing interaction with 
underdoping. Moreover, a comparison of our data to pen- 
etration depth measurements indicates that the slope of 
the superfluid density is renormalized by a doping de- 
pendent factor, implying that a non-interacting picture 
of quasiparticle excitations around the nodes of the d- 
wave order parameter is inappropriate. This has obvious 
implications for other low temperature measurements in 
the high temperature cuprate superconductors, such as 
specific heat, NMR, and microwave and thermal conduc- 
tivity, which are usually quantified by theories which do 
not take into account these renormalizations. 
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